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,_5 ' Abstract. A scheme for fast, compact, and controllable acceleration of heavy 

t^. particles in vacuum has been recently proposed [Peano F et al. 2008 New J. Phys. 

Y^ \ 10 033028], wherein two counterpropagating laser beams with variable frequencies 

• i-H ' drive a beat-wave structure with variable phase velocity, leading to particle trapping 

^,^. and acceleration. The technique allows for fine control over the energy distribution and 

J;^ I the total charge of the accelerated beam, to be obtained via tuning of the frequency 

^-f- variation. Here, the theoretical bases of the acceleration scheme are described, and the 

possibility of applications to ultrafast muon acceleration and to the prompt extraction 

of cold-muon beams is discussed. 
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1. Introduction 

The recent developments in the ultraintense-Iaser technology in the infrared (IR) optical 
regime made compact acceleration of charged particles a central research topic. The 
main research directions in the field are electron acceleration in laser-driven plasma 
waves [1-5] and proton or light-ion acceleration in laser-solid interactions (cf. [6-10] 
and, for a review, [11-13]). Typically, these acceleration schemes are indirect, in the 
sense that the acceleration is provided by laser-driven field structures in plasmas. 
As an alternative, direct, all-optical acceleration by the electromagnetic (EM) field 
of ultraintense laser pulses have been considered for electrons [14-26] and, very 
recently, also for ions [27-30]. In [29], it has been shown that the beating of two 
counterpropagating EM waves can trap and accelerate heavy particles, such as protons, 
heavy ions, and muons, even with nonrelativistic radiation intensities (e.g., -C 10^^ 
W/cm^ for proton acceleration with IR lasers), provided that at least one of the EM 
waves has variable frequency. The method offers the possibility of efficient control over 
the relevant beam features (mean energy, energy spread, and total accelerated charge), 
and can be adapted to different sources of charged particles (including beams, and 
tenuous gases and plasmas), serving either as a pre-accelerator, creating an energetic 
beam from heavy particles with low energy, or as an energy booster, accelerating 
particles from a beam with a given energy. If the necessary technological requirements 
are met [29,30], the technique can represent an alternative for compact production of 
high-quality ion beams, as well as a unique option for ultrafast acceleration of muons 
(which could be important to minimize decay losses) during, or immediately after, the 
muon cooling process [31-39] in a muon collider or a neutrino-factory machine [40-47]). 
In this Paper, the theoretical bases of the all-optical acceleration scheme are 
described (Section [2]), providing estimates for the beam properties and discussing the 
laser requirements (Section [3]), and the possibilities of applications to ultrafast muon 
acceleration and to the prompt extraction of cold muon beams are investigated (Section 

iD. 

2. Physical mechanism 

The acceleration technique presented in [29] exploits two counterpropagating EM 
waves with variable frequency in order to trap particles in a beat-wave structure with 
controllable phase velocity, according to the configuration sketched in Fig. [H The 
physical mechanism of acceleration can be analyzed in the framework of a relativistic, 
one-dimensional, single-particle theory, in which the amplitude of both lasers is constant 
(two-dimensional particle-in-cell simulations, accounting for spatial charge distributions 
and finite-size laser pulses, confirmed the validity of the present theory [29,30]). 
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Figure 1. (Color online) Sketch of the acceleration technique: a test particle is 
trapped and accelerated by the slow ponderomotive beat wave generated by two 
counterpropagating, variable-frequency lasers. The thin, solid lines represent the vector 
potential of the two lasers; the thin, dashed line represents the sum of the vector 
potentials in the superposition region; the thick, solid line represents the ponderomotive 
beat-wave. 



2.1. Field configuration and basic equations 

In the present theory, the lasers are described as two coherent EM waves (labeled with 
i = 1,2) propagating along the x direction, with vector potentials 

(1) 



Aj (x, t) = Aj< cos (6) sin [$j {^j)] ey + sin (6) cos [^j {^j)] e, 

where By and e^ are unit vectors, 6 determines the polarization type (e.g., 6 = 7i/2 for 
linear polarization in the z direction, 6 = 7r/4 for circular polarization), ^i = x — ct (laser 
1 propagates from left to right), ^2 = —x — ct (laser 2 propagates from right to left), Aj 
is the amplitude, and $j (^j) are two arbitrary functions with first derivative $j- (^j) > 
(henceforth, for the purpose of readability, the arguments ^j will be omitted unless 
necessary). The local wavenumbers, kj (x, t) = -^, and frequencies, Uj {x,t) = —-gf, 
are ki{x,t) = ^[, ui{x,t) = c^[, k2{x,t) = —^2, and LJ2{x,t) = c$2. The 
superposition of the two EM waves can be described in terms of the two corresponding 
beat waves: (i) a fast, superluminal beat wave, having phase ^F{x,t) = $1 -|- $2, 



wavenumber K{x,t) = ^ = |($'i - $2), frequency fi(x,t) = -^ = f($'i + ^ 



■27) 



and phase velocity V^{x,t) = ^ = c ^j_^f , and (ii) a slow, subluminal beat wave, 
having phase $g(x,t) = $1 — $2 wavenumber k(x,t) = -^ = |($i + $2)5 frequency 
uj{x,t) = — ^ = |($i — $2) ^^"i phase velocity v^{x,t) = f = c ^j'J . The ratio 
between the frequencies of the two laser beams in x = at t = is assumed to be such 
that the phase velocity of the slow beat wave matches the initial velocity of the particle, 
c/?o, i.e., such that v^{0, 0) = c/5o, leading to ^^|^ = j^. 

By resorting to the conservation of the transverse canonical momentum, Py = 
— - (Ai -|- A2) ■ By and pz = — - (Ai -|- A2) ■ e^, the x component of the equation of 



motion for the particle is written, in cgs units, as 



dpx 
dt 



q^ 



d 



2-fMc^ dx 



(Ai + A2 



(2) 
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where q and M are the particle charge and mass, respectively, and 7 is the Lorentz 
factor, defined as M^c^7^ = M'^c'^ + p^c^ + q^ (Ai + A2) . The squared vector 
potential in Equation ([2]) is (Ai + K2Y = {Aj + Al)/2 + AiA2Cos($s) + [sin^l^') - 
cos2(^)] [2A1A2 cos ($f) + {Al + Al) cos ($3) cos (<|)f) - {Al - Al) sin ($3) sin ($f) ]/2. In 
the special case of circular polarization, sin^ {6) = cos^ (6*), (Ai + A2)^ = | (A^ + ^2) + 
A1A2 cos($s) with arbitrary Aj for both trapped and untrapped particles. For a generic 
polarization, the fast terms containing $f can be averaged out provided that the two 
following conditions are met: (i) the variation of x and p^ associated with the transverse 
quiver motion is negligible, which is verified if Aj = jj^ <^ 1 (here, Aj = ^clj, where 
m and e are the electron mass and the elementary charge, and aj = -^ is the usual 
normalized vector potential, corresponding to the peak transverse momentum of an 
electron in the laser field); (ii) along the particle trajectory Xp(t), the frequency of 
the slow beat wave is much lower than the frequency of the fast beat wave, which is 
automatically satisfied for trapped particles, because in the comoving frame the slow 
frequency stays always close to zero [30]. Thus, Equation ([2]) can be averaged over the 
fast time scale to yield the ponderomotive equation of motion 

-^ = ir^j^ cos $1 - $2 = -^^k {x, t) sm $1 - $2 , 3 

dt 2'y ax 7 

where dimensionless quantities (denoted with hatted symbols henceforth) have been 

adopted: k = ^i with k^ = A; (0,0), t = koct, x = k^x, and Px = ^ [in these 

units, the initial frequencies of the lasers are 0)1(0, 0) = 1 + /5o and 0)2(0, 0) = 1 — /?o]- 

In Equation ([3]), the averaged Lorentz factor is defined as 7^ = 1 + p^ + Al/2 + 

A2/2 + 74iy42Cos($i — $2) [48]. Equation ([3]) can be obtained from the Hamiltonian 

je{x,p,,i) = [1 + pl + Al/2 + Al/2 + A,A2Cos{^, - ^2)f\ using ^ = -^; 

accordingly, the variation of the particle energy, -^ = ^^, is given by 

d7 A1A2 d ,^ ^ , A1A2 ^ /. -\ . /:. :. N /.S 

-J = ^r^^r cos ("^1 - "^2) = ^^^ a;, t) sm ($1 - $2) . (4) 

dt 27 C# -y ^ ' 

Equation (jl]) shows that the presence of a frequency variation is a necessary condition 
to effectively transfer energy from the beating EM waves to the particle over long times. 
If the frequency of both lasers is constant (i.e., if $'^ and $2 are constants). Equation 
(jl]) reduces to 

d7 ^ ^1^2^^ ^.^ ^2 (a; _ (3^t) + </)o] , (5) 

dt 7 

where 0o = ^i (0) — $2 (0) is the initial phase of the beat wave. In this situation, 
7 cannot grow over long times (in the frame of reference where /3o = 0, J'^ is time- 
indepedent and -^ = 0), being bound to oscillate between (F ± (3qP)/{1 — /9q), where 
F = [1 + Al/2 + Al/2 + iii2cos(0o)]^/' P = {iii2[l + cos(0o)]}^/2 [30]. However, 
if at least one of the lasers has variable frequency, causing the phase velocity of the 
ponderomotive beat wave to vary, and continuous energy transfer from the beating EM 
waves to the particle is possible. 
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2.2. Resonant phase-locked solutions 

A general resonant, phase-locked solution of Equation ([3]), X {{)., is defined by 

$1 (X it) - t) - $2 (-X it) - t) = 00. (6) 

A particular resonant solution can be found by choosing the expression for $i (or $2) 
arbitrarily and then seeking the expression of $2 (or "^i) for which Equation ([3]) and 
condition ([6]) are simultaneously satisfied. Assuming, without loss of generality, that $2 
is known, imposing condition (El), and using 7qj^ = 1 + Al/2 + A2/2 + A1A2 cos (0o) and 
T"ii = 7o±T (with T being the proper time, such that ^ = 7), Equations ([3]) and (jl]) can 
be combined (cf. [30]) to obtain the evolution equation for ^1 = X — t as 

^ = /io$2Hi-2t(r„)l, (7) 

where /iq = A1A2 sin (0o) /7ox- Through integration of Equation (jTj), one can determine 
the resonant particle trajectory in the parametric form {X (th) , t (th) } and then use 
Equation (Q to determine $1. For a generic shape of $2, this process usually requires 
a numerical approach, because of the intricate dependence of t on ^. In the special 
case in which the frequency of laser 2 is constant, i.e., $2 = 1 ~ Po^ the right-hand 
side of Equation ([7]) does not depend on t and the calculation can be performed 
analytically [30]: integration of Equation ([7]) yields ^ (Ty) = y (/iTn — 2) ry, where 
/W = /io/7oih with 7o|| = Vl - /?o> and Cq = 70,1 (1 - /3o) = VI - /^o/Vl + /^o, leading 
*o ^(^tl) = 4^[co(/^^l|-2)/^r|| - 21og(l-^r||)], t(r|,) = -^[cg (/xrn - 2) /ir,, + 
21og(l-/xr|,)],p,|(r||) = ^ = [l-c2(l-/ir||f ]/[2co(l-/ir„)], and 7,1 (r,,) = ^ = 
[1 + eg (1 - /ir„)' ]/[2co (1 - /ir,,) ]. Finally, by noticing that 1 - fir, = (1 + ^iif\ X 
and t can be written as functions of .^1 and replaced in Equation ([6]), yielding 

$1 (6) = <^o + ^ log (1 + 2/io (1 + /5o) 6) • (8) 

It is useful to analyze the behavior of the resonant solution in the early stage of 
the acceleration, t <C /i~^, and in the asymptotic limit, t ^ /i~^. For t <^ /i~^, the 
above expression for i{^) gives, to second order in i, ^(i) ^ — t — i^jjit^ and the 

explicit dependence of py, 7y, and Xy on t is written as Pwit) ~ 7oii/?o + /ut + y—iiH"^, 
7ii(t) ^ 7o|| -K l3oiJ,i + +'^011 ° ^2£2^ g^j^j ^^£^ ^ /3o^ + Y^fJ'i'^- In the nonrelativistic 
limit (/5o ^ 0, 7oy ^ 1), the particle accelerates uniformly and its energy grows 
quadratically in time as ^filt"^- For t ^ fi"^, the behavior of the solution depends 
on whether laser 1 and the particle are copropagating (/i > 0) or counterpropagating 
(yu < 0): if /i > 0, T\\(t) ^ -[1 — exp (^ — 2co/it)] and, consequently, X(t) ^ t — |2- and 

7y(t) ~ Pwii) ~ 2c~ ^-^P (^'^0'^^ ~" ^)' if /^ < 0, r||(t) :^ 2(t/|/iCo|) and, consequently, 
X(t) ^ — t and 7||(t) ~ Pwit) ^ (|/iCo|t) . The behavior is asymmetric because, if 
/i > 0, the frequency must be increased to maintain phase- locking, causing a continuous 
increase of the ponderomotive force, whereas, if /i < 0, the frequency must be decreased, 
causing a continuous decrease of the ponderomotive force. 
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2.3. Particle trapping 

Effective particle acceleration can occur also when $i and $2 do not match the resonant 
solution exactly. In fact, it is sufficient to choose $1 and $2 appropriately in order 
to guarantee that the particle trajectory stays close to the beat-wave trajectory until 
reaching the desired energy, i.e., that the particle is trapped by the beat wave. For 
given $1 and $25 the beat-wave trajectory x^^^ (t) (namely, the trajectory of the point 
of the beat wave having constant phase 0o) is determined by solving the equation 
•^1 {X(j,Q — t) — $2 {~i<j>o ~ t) = 00 with respect to x^^q [for resonant trajectories, 
X (t) = XtjjQ (t)]. In order to determine the conditions for the occurrence of trapping. 
Equation ([3]) is expressed in terms of the phase difference between particle and wave, 
i/j = 2[x — X(^(, (t) ], and using the proper time r, as 



-U{^,t) 



(9) 



dr^ dijj 

with U iij, t) = 2iii2 cos {<l>i[ei,0o^(T) + f ] - $2[6,</-o (j) - f [} + 2«0o (^) ^, where^ei,0o. 
^2,</,o, and a^g depend on r through i as ^i<^„ (r) = x^^ (i) -t, ^2,9^0 (^) = "^^o (^) -i^ and 
Oi<f,o (r) = ^x^o (t) = 7 (t) f^ [7 (t) f^x^g (t)]. According to Equation ([9]), trapping is 
allowed only if the frequency variation is slow enough to guarantee that the effective 
potential U presents local minima. For trapped particles, IV'I <C '^^j<j)o and U can be 




T\v^® 



Figure 2. (Color online) Space-time evolution of the dimensionless ponderomotive 



force 



dx 



COS ($1 — $2) (color scale, red and blue denoting accelerating and 



decelerating regions, respectively), and trajectories of a resonant particle (dashed 
line), a non-resonant, trapped particle (green spheres), and a non-resonant, untrapped 
particle (yellow spheres), for a situation in which /?o = 0, $2 = 1, and $1 is given by 
expression ^, with A1A2 = 5 x 10~^ and ^0 = ^tt/G. The inset represents the effective 
potential U at t — 0, gray marking the trapping region, and the corresponding initial 
positions of the three particles. 



approximated as U (-0, r) ~ 2A1A2 cos {k[x^Q (i) , i\ijj + 0o} + Sa^^ (t) ^p leading to the 
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necessary condition for particle trapping 

1 000 (t) I < Om [i] = AiA2k [x^o (t) , i] , (10) 

where Om [t] is the maximum value of the ponderomotive force. Thus, trapping is 
allowed only if the maximum ponderomotive force of the beat wave is greater than 
the inertial force associated with the beat-wave acceleration. The particular situation 
in which |a;0„ (t) | = Om (i) corresponds to a resonant solution in the particular case 

00 = ±vr/2. When $i and $2 correspond to a resonant solution and |0o| 7^ ^/2, regions 
where trapping is possible always exist, although the solution is stable for cos (0o) < 
and unstable for cos (0o) > 0, with ip = corresponding to the bottom and top of a 
potential well, respectively. The relative width of the potential well in t/(^,0), Aip, 
can be used to estimate the trapping efficiency rj (defined as the fraction of trapped 
particles) as r^ ^ ^ ~ 1 — ^ arcsm[a^g{0) /aM]- Examples of typical trapped and 
untrapped trajectories are shown in Fig. [2J 

2.4- Acceleration with linearly chirped lasers 

An important case of variable-frequency lasers, particularly relevant in experiments, 
is that of linearly chirped beams, in which $1 and $2 take the form $1(^1) = 

001 + (1 + /3o)6 + (^i^i and $2(6) = 0O2 + (1 - /^o)^ + 0-2^2^ where 0oj are constant 
phases and ctj are the chirp coefficients. The local frequencies and wavenumbers are 
given by u)i(x, i) = ki{x, i) = 1 + /3o + 2ai^i and 0)2(3;, i) = -k2{x, i) = 1 - /3o + 2cr2^2- 
The wavenumber and the frequency of the slow beat wave are k{x,t) = 1 -|- a_x — cr+t 
and uj{x,t) = Po + a+x — cr_t, where (T_ = cti — (J2 and a^ = di -\- 02- The phase of 
the ponderomotive beat wave is $1 — $2 = 0o + 2(x — %t — o^xt) + cr_(a;^ + t^), with 
00 = 0oi~0o2- Consequently, the beat- wave trajectory is X(^o(t) = ^^t — ^{1— [1— 2(cr+ — 
/3ocr_)t + (a^ — cr^)t^] }. If only one laser is chirped, i.e., if (72 = or a\ = 0, one finds 

£^l^{i)=i-j-^{l-[l-2a,il-Po)tr^'}andx^^^^{i) = -Uj;{l-[l-2a2il + m'^"}, 
respectively. This suggests that the acceleration can be improved by chirping both lasers, 

with o"iO"2 < 0. However, when the particle and the chirped laser are copropagating (e.g., 

cTi < and (T2 = 0), a^^ (t) is monotonically decreasing in time, causing the trapping 

regions to widen [30]; on the contrary, when the particle and the chirped laser are 

counterpropagating (e.g., o"i = and 0-2 > 0), a^f.^ (i) is monotonically increasing in 

time, causing the trapping regions to narrow. When chirping both lasers (cti 7^ and 

c"2 7^ 0), the use of a stronger chirp in the laser copropagating with the particle is in 

general convenient: in such a situation, the general expression for a^^ (t) has a minimum 

for a given time (which is t = if 0"i = — cr2), such that the width of the trapping regions 

increases until ctt^g (t) reaches its minimum, and starts decreasing immediately after. 

This is illustrated in Fig. [3l showing the trapping and acceleration process for a case 

wherein only one laser is chirped (case A) and for a case wherein both lasers are chirped 

(case B): in case A, the amplitude of the trapping region, A-0, grows monotonically 

towards its maximum value, Atp = 27r, and the acceleration continues indefinitely; in 
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case B, A-?/; reaches a maximum and then starts decreasing until synchronization with 
the beat wave is lost and the acceleration stops. 




5 10 

Time [ps] 



Figure 3. (Left panel) Width of trapping region as a function of time and (right 
panel) evolution of particle energy as a function of the acceleration distance when only 
laser one is chirped with ai = —A1A2 (case A), and when both lasers are chirped with 
CTi — —0.65^1^2 and (T2 = 0.35^1^2 (case B). Particle type: muon. Initial particle 
velocity: /3o = 0. Initial laser intensities: Ii = I2 = 1-3 x IQ-'^^W/cm'^. Initial laser 
wavelenghts: 800 nm. 



3. Beam properties and laser requirements 

Basic estimates for the principal acceleration features are readily obtained in the 
nonrelativistic regime, wherein $j can always be approximated as ^j{^j) ~ 0oj + ^j + 
crjC^i with linear-chirp coefficients aj = \^'j{^j)- To leading order in t, this leads to 
the beat-wave trajectory x^g{t) ~ |cr_t^. Hence, the momentum and kinetic energy 
of trapped particles can be approximately expressed as a_t and |cr^t^, respectively. 
The maximum energy gain, AeM, which occurs when using the maximum chirp that 
allows for trapping, |cr_| = A1A2, can then be estimated as AeM[MeV] ~ 0.8 g^[e] 
M-3[amu] /i[1020W/cm2] JailO^OW/cm^] A^[/xm] AT2[ps], where Aq = 27r/fc(0,0) is the 
initial laser wavelength; the corresponding acceleration distance scales as Aa;[/im] f» 6 



.1/2 



M-2[amu] /i^^[1020W/cm2] I^^^[10'^°W/cm^] Ao[/im] AT2[ps]. According to these 



estimates, if sufficiently intense lasers are employed, the method can provide ultrafast 
acceleration of heavy particles over short distances. 

The maximum energy gain can also be expressed in terms of the laser-pulse energies 
as A^,,[MeV] ^ 0.08 q^[e\M-^[amn\£i[J]£2[3]\l[fxm]W^'^[fxm]. Hence, if Wq and Aq are 
given, ASm depends only on the energy of the two laser pulses, and the combination 
of pulse lengths and intensities can be tuned in order to obtain the desired balance 
between the relative width of the final energy spectrum (which scales as ^^ ^ ) and the 
total accelerated charge (which scales as ^1^42 if the laser durations are longer than the 
diffraction length, being constant otherwise). Provided that the available laser energy is 
sufficient, the other fundamental requirement of the acceleration method is the frequency 
bandwidth. The frequency excursions, AcUj, necessary to accelerate a particle from the 



:iii 
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initial velocity (3q to the final velocity (3 are determined by 
2/?o + /\ui - Auj2 
2 + AuJi + Auj2 

When only one laser is chirped, Auj are given by Ac^i = 2 {P — Pq) / {1 — P) and 
Au!2 = 2{P — Po) / {1 + P). The minimum frequency excursion per laser is obtained 
when Acui = —Auj2 = P — Po- 

The present acceleration technique remains effective also in more general situations 
involving distributed particle sources and finite-size laser pulses, provided that the 
particle density, n, is sufficiently low to prevent significant distortions of the EM beat- 
wave structure by the space-charge field, and the longitudinal and transverse variations 
in the radiation intensity are sufficiently smooth to guarantee the validity of the plane 
wave approximation in the focal region (as confirmed by two-dimensional particle-in- 
cell simulations [29,30], performed with the Osiris 2.0 simulation framework [49]). 
These considerations lead to the constraint on the particle density n[10^^cm~^] <^ 
q~^[e]Ii [10^°W/cm^]/2 [10^°W/cm^], and, to the estimate for the maximum trapped 
charge (5m[pC] ~ 16r7g[e] n[10^^cm~^] W(f[/xm]/Ao[Aim], where Wq is the overlapped spot 
size of the lasers and rj is the fraction of trapped particles. 



4. Ultrafast muon acceleration 




5 10 15 

Time [ps] 



10 15 
Time [ps] 



Figure 4. Example of single-stage muon acceleration: (a) laser frequencies along the 
muon trajectory, (b) muon velocity and beat-wave phase velocity, (c) muon trajectory, 
and (d) evolution of the muon energy. Laser parameters: Ii = I2 = lAx 10^^ W/cm^, 
CTi = —4.4 X 10^^, CTi == 4.4 X 10^^. Initial laser wavelenghts: 800 nm. 



Gaining the ability to produce high-quality muon beams would provide 
unprecedented capabilities for fundamental-physics experiments in pioneering research 
topics, ranging from investigations on neutrino-oscillation physics in neutrino-factory 
machines to energy- frontier studies in TeV-level muon-antimuon colliders [40-43]. 
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From the point of view of technology, muon acceleration is an extremely challenging 
task, involving a series of delicate processes, namely, producing high-energy muons 
via interaction of intense proton beams with appropriate targets, cooling the muon 
population in order to lower the emittance to acceptable levels (proposed methods are 
ionization cooling [31-36] and frictional coohng [37-39]), and, finally, accelerating cold 
muons to high energy. Moreover, due to the short muon lifetime (2.2 /is), these processes 
must be performed rapidly to minimize decay losses. 

The all-optical acceleration scheme described here, as well as other conceivable 
laser-based techniques, might be appealing because of its potential for ultrafast 
acceleration, which could be particularly important in the first stages of muon 
acceleration, especially during, or immediately following, the cooling process. In the 
example shown in Fig. HI a test muon is accelerated from zero velocity to 0.9c in just 25 
ps over 4 mm in a single stage. Unfortunately, such an acceleration pattern would require 
very large frequency excursions, which are extremely difficult to obtain, at least with 
current laser technology. A possible strategy to circumvent this obstacle is dividing the 
acceleration process in several stages such that the acceleration direction in each stage 
is always perpendicular to the propagation direction of the beam at the entrance of the 
stage [50]. Theoretically, this would permit the achievement of arbitrary large A/3, by 
repeatedly using a limited Acu. An example of such a multi-stage acceleration pattern 
is depicted in Fig. [5] where moderate-intensity lasers and 10 stages with frequency 
variations on the order of 50 % are used to accelerate a muon from zero- velocity to 0.7c. 
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Figure 5. Multi-stage muon acceleration: (left panel) muon trajectory in the x- 
y plane and (right panel) evolution of the muon energy (solid line) during multi- 
stage acceleration. The duration of each stage is for 4 ps. In each stage s, the laser 
propagation direction is perpendicular to the muon momentum at the entrance of the 
stage, pos (the values of the angle between pos and the e^: are indicated) . Laser 1 is 
linearly chirped with ais — 5.47 x 10^^/^/1+^08' while the frequency of laser 2 is 
fixed; the laser intensities are /i = /2 = 5 x 10^^ W/cm^. In the right panel, the dashed 
line refers to the evolution of muon energy corresponding to a single-stage process with 
(Ti = 5.47 X 10^^ and same laser intensities. 



Finally, another application of the present method might be the prompt extraction 
of low-temperature muons from a gas cell used for frictional cooling [37-39], in which 
the muon distribution relaxes to thermal equilibrium with energy below the ionization 
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threshold of the gas atoms. A possible strategy could be trapping the muons in 
the beat-wave structure generated by non- ionizing (e.g., with high power and long 
wavelength) EM waves. The assessment of the possibility of actually meeting the 
necessary combination of EM-wave parameters (i.e. wavelength, intensities, and chirp 
law) will be object of future investigation. 

5. Conclusions 

All-optical trapping and acceleration of heavy particles in vacuum can be achieved by 
employing counterpropagating, intense laser beams with variable frequencies, which also 
offers unique control capabilities on important beam features such as the final energy 
and energy spread [29,30]. Allowing for ultrafast acceleration, the technique could 
provide a route towards prompt muon acceleration in the first acceleration stages in a 
muon collider or a neutrino-factory machine. 
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